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Abstract. We consider a scattering set-up with an entangler and beam splitter where 
the current noise exhibits bunching behavior for electronic singlet states and antibunch- 
ing behavior for triplet states. We show that the entanglement of two electrons in the 
double-dot can be detected in mesoscopic transport measurements. In the cotunnel- 
ing regime the singlet and triplet states lead to phase-coherent current contributions 
of opposite signs and to Aharonov-Bohm and Berry phase oscillations in response to 
magnetic fields. We analyze the Fermi liquid effects in the transport of entangled elec- 
trons. 



1 Introduction 

The availability of pairwise entangled qubits - Einstein-Podolsky- Rosen (EPR) 
pairs [1] - is a necessary prerequisite in quantum communication [2] . The prime 
example of an EPR pair considered here is the singlet /triplet state formed by 
two electron spins [3] , [4] . Its main feature is its non- locality: If we separate the 
two electrons from each other in real space, their total spin state can still remain 
entangled. Such non-locality gives rise to striking phenomena such as violations 
of Bell inequalities and quantum teleportation and has been investigated for 
photons [5], [6], but not yet for massive particles such as electrons, let alone in 
a solid state environment. In this work we discuss specific properties of trans- 
port and noise of entangled electrons as a result of two-particle coherence and 
nonlocality. 

In Sect. (2) we propose and analyze an experimental set-up (see Fig. la) 
by which the entanglement of mobile electrons can be detected in noise mea- 
surements with a beam splitter [7]. The entangler is assumed to be a device by 
which we can generate entangled electron states, a specific realization being the 
double-dot system [3] . The presence of a beam splitter ensures that the electrons 
leaving the entangler have a finite amplitude to be interchanged. Thus we can 
expect that the current-current correlations (noise) measured in leads 3 and/or 
4 are sensitive to the symmetry of the orbital part of the wave function [8] . Since 
the spin singlet of two electrons is uniquely associated with a symmetric orbital 
wave-function, and the triplet with an antisymmetric one we have a means to 
distinguish singlets from triplets through a bunching or antibunching signature. 
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Fig. 1. a) The setup for measuring noise correlations of entangled states. Uncorrelated 
electrons are fed into the entangler through the Fermi leads 1' and 2'. The entangler is 
a device that produces pairs of entangled electrons and injects one of the electrons into 
lead 1 and the other into lead 2. The entanglement can then be detected by performing 
an interference experiment using a beam splitter, b) Double-dot (DD) system contain- 
ing two electrons and being weakly coupled to metallic leads 1,...,4, each of which being 
at the chemical potential fii,...,fi4. The tunneling amplitudes between dots and leads 
are denoted by F, F. The tunneling (t) between the dots results in a singlet (triplet) 
ground state. The closed tunneling path between dots and leads 1 and 2 encloses the 
area A. 

It is well-known [10] that bosons (fermions) show bunching (antibunching) be- 
havior [11]. Antibunching is so far considered for electrons in the normal state 
both in theory [12], [13] and in experiments [14]. However, this classical effect is 
independent of phase coherence [15] and should be carefully distinguished from 
the two-particle phase-coherent effect which we propose here. 

The scheme we propose in Sect. (3) [16] consists of two coupled quantum 
dots (DD) which themselves are weakly coupled in parallel to two leads 1 and 
2 (see Fig. lb). This results in a closed loop, and applying a magnetic field, an 
Aharonov-Bohm (AB) phase ip will be accumulated by an electron traversing 
the DD. In the Coulomb blockade (CB) regime we find that due to cotunneling 
[17] the current depends on the state of the DD: the AB oscillations for singlet 
and triplets have opposite sign. The amplitude of the AB oscillations provides 
a measure of the phase coherence of the entangled state, while the period - via 
the enclosed area A - provides a measure of the non-locality of the EPR pairs. 
The triplets themselves can be further distinguished by applying a directionally 
inhomogeneous magnetic field which adds a Berry phase [18] leading to beating. 

Finally, in Sect. (4) we would like to address the following question [9]: Is 
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it possible to use mobile electrons, prepared in a entangled spin state, for the 
purpose of quantum communication? Without spin-dependent interaction we 
know that the total spin must be conserved even if the two electrons strongly 
interact with the other electrons in the mcsoscopic environment (and among 
themselves) via Coulomb interaction. It is thus not unreasonable to expect that 
we still find some spin correlations between initial and final states. But how 
much is it? And why and how do we loose some of the correlations? 

2 Noise of Entangled Electrons: Beam Splitter Set-up 

Below, we extend the standard scattering matrix approach [12] to a situation 
with entanglement. We start by writing the operator for the current carried by 
electrons with spin cr in lead a of a multiterminal conductor as 

w = E E 4a(^)«7a(^')^37 - . (1) 

j4^^ = SctjsSctj — Sf^pSctj ■ (2) 

where al^^{e) creates an incoming electron in lead a with spin a and energy e, 
and wc assume that the scattering matrix Safj is spin- and energy- independent. 
Note that since we are dealing with discrete energy states here, we normalize 
the operators aa{e) such that [aaai^), «/3(t'(^')^] = Saa'^a0Se,e' /i^, where the Kro- 
ncckcr symbol S^^^i equals 1 if e = e' and otherwise. Therefore wc also have 
included the factor in the definition of the current, where u stands for the 
density of states in the leads. We will also assume that each lead consists of only 
a single quantum channel; the generalization to leads with several channels is 
straightforward but is not needed here. 

We restrict ourselves here to unpolarized currents, 1^ = laa- The spectral 
density current fluctuations (noise) Sla = — between the leads a and /? 
are defined as 

Sa0{u>)= hm ^ / dtexp{iwt){^\SIa{t)Sl0{O)\^) , (3) 

where \^) is the quantum state of the system to be specified next^. We will now 
investigate the noise for scattering with the entangled incident state 

|±) = (4^(e2)«I^(ei) ± 4^(e2)at^(ei)) |0) . (4) 

The state |— ) is the spin singlet, \S), while |+) denotes the Sz = triplet |To)^. 
Substituting |±) defined in (4) for and using the fact that the unpolar- 
ized currents are invariant when all spins are reversed, the expectation value 

^ Note that since \'P) in general does not describe a Fermi liquid state, it is not possible 

to apply Wick's theorem. 
^ All three triplets, \To), ITT); I ii); give the same result. 



4 



Eugene V. Sukhorukov, Daniel Loss, and Guido Burkard 



{±\6IaSI^\±) can be expressed as the sum of a direct and an exchange term, 

{±\5Io,Sl0\±) = {U \SIaSIp\ U) ± (Ti \SIaSIp\ it) , (5) 

where the upper (lower) sign of the exchange term refers to triplet (singlet). 
Direct calculation of (5) gives the following result for the zero-frequency (w = 0) 
correlations, 



7(5 



A'^.A^ m /) ( A°' A^ A°' A^ \\ 

^75^(57 + "ei,e2 I ^12^21 + ^21^12/' j 



(6) 



where Y^'^^ denotes the sum over 7 = 1,2 and all 5 

We apply formula (6) now to the set-up shown in Fig. la involving four leads, 
described by the single-particle scattering matrix elements, Szi = S42 = r, and 
S41 = S32 = t, where r and t denote the reflection and transmission amplitudes 
at the beam splitter, respectively. We assume that there is no backscattering, 
S12 = S34 = Saa = 0. The unitarity of the s-matrix implies |rp + |tp = 1, and 
Re[r*f] = 0. Using (2) and (6), we obtain the final result for the noise correlations 
for the incident states |±), 

Ss3 = SM = -Su = eF\{I)\ , 

F = 2eT(l-T)(lT<5ei,.J , (7) 

where |(/)| = e/hu is the average current in all leads, T = \t\'^ is the probability 
for transmission through the beam splitter, and F is the Fano factor. Note that 
the total current 6I3 + 6I4 does note fluctuate, i.e. -|- S44 + 26*34 = 0, since 
we have excluded backscattering. 

Above results (7) imply that if two electrons with the same energies, ei =62, 
in the singlet state \S) = |— ) are injected into the leads 1 and 2, then the zero 
frequency noise is enhanced by a factor of two, F = 4eT(l — T), compared to the 
shot noise of uncorrclatcd particles, F = 2eT{l — T). This enhancement of noise 
is due to bunching of electrons in the outgoing leads, caused by the symmetric 
orbital wavefunction of the spin singlet |5). On the other hand, all three triplets 
1+) exhibit an antihunching effect, leading to a complete suppression of the 
zero-frequency noise, Saa = 0. The noise enhancement for the singlet \S) is a 
unique signature for entanglement (there exists no unentangled state with the 
same symmetry) , therefore entanglement can be observed by measuring the noise 
power of a mesoscopic conductor as shown in Fig. la.^ 



3 Probing Entanglement of Electrons in a Double-Dot 

The DD system (see Fig. lb) contains 4 metallic leads which arc in equilibrium 
with associated reservoirs kept at the chemical potentials jii, i — 1, ... ,4. The 

^ These results remain valid in the presence of a Fermi sea. 



Transport and Noise of Entangled Electrons 



5 



leads are weakly coupled to the dots with tunneling amplitudes F and F, and 
the leads 1,2 are coupled to both dots and play the role of probes where the 
currents are measured. The leads 3 and 4 are feeding electrodes to manipulate 
the electron filling in the dots. The quantum dots contain one (excess) electron 
each, and are coupled to each other by the tunneling amplitude t, which leads 
to a level splitting [3], [4] J = Et - At^ /U in the DD, with U being the 

single-dot Coulomb repulsion energy, and E^/^ are the singlet /triplet energies. 
We recall that for two electrons in the DD (and for weak magnetic fields) the 
ground state is given by a spin singlet. For convenience we count the chemical 
potentials fii from Eg. The coupling F to the feeding leads can be switched off 
while probing the DD with a current. Here we assume that F = Q. 

Using a standard tunneling Hamiltonian approach [19], we write H = Hq + V , 
where the first term in _ffo = -^d + H1+H2 describes the DD and Hi 2 the leads 
(assumed to be Fermi liquids). The tunneling between leads and dots is described 
by the perturbation V = Vi + V2, where 



K = r ^ [Dl^Cn,s + clsDn,s\ , Dn,s = e±^^/^rfl,« + e^^^/^ds.a , (8) 



and where the operators Cn,s and d„,s, n = 1,2, annihilate electrons with spin 
s in the nth lead and in the nth dot, resp. The Peierls phase in the hopping 
amplitude accounts for an AB or Berry phase (see below) in the presence of a 
magnetic field. The upper sign belongs to lead 1 and the lower to lead 2. Finally, 
we assume that spin is conserved in the tunneling process. For the outgoing 
currents we have /„ = ieFj^s [-^n.sCn.s — cJi.s^'i s] • The observable of interest 
is the average current through the DD system, / = (/2). 

From now on we concentrate on the CB regime where we can neglect double 
(or higher) occupancy in each dot for all transitions including virtual ones, i.e. 
we require /xi,2 < U. Further we assume that /xi,2 > J, ksT to avoid resonances 
which might change the DD state. The lead-dot coupling F is assumed to be 
weak so that the state of the DD is not perturbed; this will allow us to retain only 
the first non- vanishing contribution in F to /. Formally, we require J > 2TTi'tF^, 
where vt is the tunneling density of states of the leads. In analogy to the single- 
dot case [17] we refer to above CB regime as cotunneling regime. 

Continuing with our derivation of /, we note that the average (...) = Trp {. . .} 
is taken with respect to the equilibrium state of the entire system set up in the 
distant past before V is switched on [19]. Then, in the interaction picture, the 
current is given by 



The leading contribution in F to the cotunneling current involves the tunneling 
of one electron from the DD to, say, lead 2 and of a second electron from lead 1 
to the DD (see Fig. lb). This contribution is of order V2V1, and thus I oc F'^, as 



s 




(9) 
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is typical for cotunneling [17]. Taking the trace over Fermi leads, we arrive then 
at the following compact expression for the cotunneling current 

A,j = fii~ ^12 + E,- Ef . (10) 

This equation shows that in the cotunneling regime the initial state \i) (with 
weight Pi) of the DD is changed into a final state |/) by the traversing elec- 
tron. However, due to the weak coupling F, the DD will have returned to its 
equilibrium state before the next electron passes through it. 

For small bias, — ^2! < J, only elastic cotunneling is allowed, i.e. Ei = Ef. 
However, this regime is not of interest here since singlet and triplet contributions 
turn out to be identical and thus indistinguishable. We thus focus on the opposite 
regime, — /X2I > J, where inelastic cotunneling* occurs with singlet and 
triplet contributions being different. In this regime we can neglect the dynamics 
generated by J compared to the one generated by the bias ("slow spins"), and 
drop the energies Ei and Ef in (10). Finally, using 1 = |/)(/| we obtain 

I = eni.^r^C{^)^^^^ , (11) 

^du') + cosip{d\^,dudlgd2s')\ ■ (12) 

s,s' 

For the purpose of our analysis we assume that the DD is in its ground state. 
Equation (12) shows that the cotunneling current depends on the properties of 
the ground state of the DD through the coherence factor C{(p) given in (12). 
The first term in C is the contribution from the topologically trivial tunneling 
path which runs from lead 1 through, say, dot 1 to lead 2 and back. The second 
term (phase-coherent part) in C is the ground state amplitude of the exchange 
of electron 1 with electron 2 via the leads 1 and 2 such that a closed loop is 
formed enclosing an area A (see Fig. lb). Thus, in the presence of a magnetic 
field B, an AB phase factor cp = ABe/h is acquired. 

Next, we evaluate C{ip) explicitly in the singlet-triplet basis. Note that only 
the singlet \S) and the triplet \To) (see (4)) are entangled EPR pairs while 
the remaining triplets |T_|_) = | ||), and |T_) = | J,J,) are not (they factorize). 
Assuming that the DD is in one of these states we obtain the important result 

^ f 2 - COS , for singlet 

12 + cos If , for all triplets . ' 

Thus, we see that the singlet and the triplets contribute with opposite sign to 
the phase-coherent part of the current. One has to distinguish, however, carefully 

"* Note that the AB effect is not suppressed by this inelastic cotunneling, since the 
entire cotunneling process involving also leads is elastic: the initial and final states 
of the entire system have the same energy. 
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the entangled from the non-entangled states. The phase-coherent part of the en- 
tangled states is a genuine two-particle effect, while the one of the product states 
cannot be distinguished from a phase-coherent single-particle effect. Indeed, this 
follows from the observation that the phasc-cohcrcnt part in C factorizcs for the 
product states T± while it does not so for S, Tq. Also, for states such as | ti) the 
coherent part of C vanishes, showing that two different (and fixed) spin states 
cannot lead to a phasc-cohcrcnt contribution since wc know which electron goes 
which part of the loop. Finally we note that due to the AB phase the role of the 
singlet and triplets can be interchanged which is to say that we can continually 
transmutate the statistics of the entangled pairs S, Tq from fcrmionic to bosonic 
(like in anyons): the symmetric orbital wave function of the singlet S goes into 
an antisymmetric one at half a flux quantum, and vice versa for the triplet Tq. 

We would like to stress that the amplitude of the AB oscillations is a direct 
measure of the phase coherence of the entanglement, while the period via the 
enclosed area A = h/eBo gives a direct measure of the non-locality of the EPR 
pairs, with Bq being the field at which = 1. The triplets themselves can be 
further distinguished by applying a directionally inhomogeneous magnetic field 
(around the loop) producing a Berry phase [18], which is positive (negative) 
for the triplet m = 1(— 1), while it vanishes for the EPR pairs 5*, Tq. Thus, we 
will eventually see beating in the AB oscillations due to the positive (negative) 
shift of the AB phase ^ by the Berry phase, (p = $±^^. 

4 Transport of Entangled Electrons 

We consider the general scenario of the transport of entangled electrons in a 
mesoscopic system [9]. In a first step we inject entangled electrons into the 
leads and create the state |12) = I±) (see (4)) on the top of the Fermi sea (as 
discussed e.g. in Sect. (2), see Fig. la). In a second step, we perform a quantum 
measurement of the state. As a measure of correlations we consider transition 
amplitudes between an initial and a final state. We begin with the simplest case 
given by the wave function overlap of |12) with |34), 



where the upper (lower) sign refers to triplet (singlet). If the quantum numbers 
coincide, 1 = 3, and 2 = 4, the overlap assumes its maximum value 1, reflecting 
maximum correlation between the two states. 

Next we generalize this concept to leads which contain many interacting 
electrons besides the two entangled electrons. We use a similar overlap as a 
measure of how much weight remains in the final state |34, t) when we start from 
some given initial state |12). The overlap (14) now becomes a triplet-triplet or 
singlet-singlet correlation function 



(12|34) = S13S24 T '^i4'^23 , 



(14) 



G*/^(12, 34; t) = -G(l, t) G(2, t) [6^3624 T S14S23] 



(15) 
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where we have assumed that there is no interactions between lead 1 and 2. Thus 
the problem is reduced to the evaluation of (time-ordered) single-particle Green's 
functions G{l,t), G{2,t) pertaining to lead 1 and 2, resp. (these leads are still 
interacting many-body systems though). 

For the special case t = 0, and no interactions, we have G = — i, and thus 
G^/^ reduces to the rhs of (14). For the general case, we evaluate G close to the 
Fermi surface and get the standard result [20] 

G{e, t) « -iz,0{e - eF)e-^^*-^'* , (16) 

where e is the quasiparticle energy, ep is the Fermi energy, and l/A is the 
quasiparticle lifetime. In a 2DEG, oc (e — ep)^ log(e — ep) [21] within the 
random phase approximation (RPA). Thus, the lifetime becomes infinite when 
the energy of the added electron approaches ep. 

Now, we come to the most important quantity in the present context, the 
quasiparticle weight, z-p = z^^, evaluated at the Fermi surface; it is defined by 

(17) 

where S{e, u)) is the irreducible self-energy occurring in the Dyson equation. The 
quasiparticle weight, < < 1, describes the weight of the bare electron in the 
quasiparticle state e, i.e. when we add an electron with energy e > ep to the 
system, some weight (given by 1 — z^) of the original state e will be distributed 
among all the electrons due to the Coulomb interaction. 

Restricting ourselves now to energies close to the Fermi surface we have 

G*/^(12,34;t) = z|[(5i3(524T<5i4M , (18) 

for all times satisfying < t ^ l/T^- Thus we see that it is the quasiparticle 
weight squared, Zp, which is the measure of our spin correlation function G*/*^ 
we were looking for. It is thus interesting to evaluate zp explicitly. This is indeed 
possible, again within RPA, and we find after careful calculation [22] 

zp = l- r.(i + -) , (19) 

Z TT 

in leading order of the interaction parameter Ts = l/gpae, where ae = e^h^ jme? 
is the Bohr radius. In particular, in a GaAs 2DEG we have ob — 10.3 nm, 
and = 0.614, and thus we obtain from (19) the value zp = 0.665. We note 
that a more accurate numerical evaluation of the exact RPA self-energy yields^ 
Zp = 0.691 [22]. Thus, we sec that the spin correlation is reduced by a factor of 
about two as soon as we inject the two electrons into separate leads consisting 
of interacting Fermi liquids in their ground state. 

^ For 3D metallic leads with, say, rs = 2 (e.g. r^" = 2.67) the loss of correlation is 
somewhat less strong, since then the quasiparticle weight becomes zf = 0.77 [19]. 



Zf 



1 - —ReS(eF,LO = 0) 
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